IPMU09-0066 



Dark matter as integration constant in 
Hofava-Lifshitz gravity 

Shinji Mukohyama 

Institute for the Physics and Mathematics of the Universe (IPMU) 

The University of Tokyo 
5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8582, Japan 

Abstract 

In the non-relativistic theory of gravitation recently proposed by Hofava, 
the Hamiltonian constraint is not a local equation satisfied at each spatial 
point but an equation integrated over a whole space. The global Hamilto- 
nian constraint is less restrictive than its local version, and allows a richer set 
of solutions than in general relativity. We show that a component which be- 
haves like pressureless dust emerges as an "integration constant" of dynam- 
ical equations and momentum constraint equations. Consequently, classical 
solutions to the infrared limit of Horava-Lifshitz gravity can mimic general 
relativity plus cold dark matter. 



1 Introduction 



Dark energy and dark matter are two major mysteries in modern cosmology. As- 
suming that general relativity is correct at long distances up to cosmological scales, 
precision observational data indicates that more than 90% of our universe consists 
of dark energy and dark matter. Although some gravitational properties of the dark 
components are known, they are not optically observed and, thus, we do not know 
what they really are. This situation makes us suspect that modifying gravity in the 
infrared (IR) might address the mysteries of dark energy and/or dark matter. 

Recently a power-counting renormalizable theory of gravitation was proposed 
by Horava [H [2]. One of the most important aspects of the theory is that in the 
ultraviolet (UV) it is fundamentally non-relativistic and exhibits the Lifshitz scale 
invariance 



with dynamical critical exponent z = 3. Horava's theory is considered as a poten- 
tial candidate for the theory of quantum gravity and is often called Hofava-Lifshitz 
gravity. Various aspects of this theory have been investigated [5]-[in]. 

Hofava-Lifshitz gravity has not yet been intended to be a unified theory. Clearly, 
further developments or/and embedding into a "bigger" theory is needed. For ex- 
ample, since the "limit of speed" is an emergent quantity in the IR, different species 
including those in the standard model of particle physics must be related to each 
other in the framework of Horava's theory so that the "limits of speed" for different 
species in the IR agree with the "velocity of light" B This obviously indicates that 
embedding of this theory into a unified theory (or other way around) is necessary for 
the theory to be a part of the real world. 

Still, it is interesting to investigate universal properties of the theory § and its 
cosmological implications, in parallel with those fundamental issues. For example, the 
z = 3 Lifshitz scaling not only is the origin of the power-counting renormalizability but 
also leads to a number of interesting cosmological consequences, such as generation of 
scale-invariant cosmological perturbations from a non-inflationary epoch of the early 
universe [7] and a particular scaling of radiation energy density (oc a~^) |iT^ . 

The purpose of the present paper is to point out that Horava-Lifshitz gravity can 

^ Note, however, that renormalizabiUty has not yet been estabhshed in a rigorous manner beyond 
the level of power-counting. 

^See e.g. refs. [41]|42lll3j for tight experimental limits on Lorentz violation. 

^In this respect, the so called detailed balance condition is neither essential nor universal for 
Hofava-Lifshitz gravity as already stated in [T]. 
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mimic general relativity plus cold dark matter. 

2 Basic idea 

Before explaining why the IR limit of Hofava-Lifshitz gravity can behave like gen- 
eral relativity plus cold dark matter, let us remind ourselves about the structure of 
Einstein's general relativity since the existence of dark matter was suspected by as- 
suming general relativity. General relativity fully respects 4-dimensional spacetime 
diffeomorphism invariance as the fundamental symmetry of the theory. As a result, 
it has four constraint equations: one called Hamiltonian constraint and three called 
momentum constraint. These constraints must be satisfied at each spatial point at 
each time. However, since the constraint equations are preserved under time evo- 
lution by dynamical equations, i.e. other components of the Einstein equation, it 
is also possible to impose the constraint equations only on an initial hypersurface 
and to solve dynamical equations afterwards. In this case, constraint equations are 
automatically satisfied at late time. 

As an illustration, let us consider a flat Friedmann-Robertson- Walker (FRW) 
spacetime driven by components with equations of state Pi = Pi{pi), where pi and 
Pi are energy density and pressure of the i-th component. Because of the spatial ho- 
mogeneity, the momentum constraint is trivial. On the other hand, the Hamiltonian 
constraint gives the famous Friedmann equation: 

. 2 n 

3^ = 87rGjvEPi' (2-1) 

where a is the scale factor of the universe, a dot represents time derivative and n is 
the number of components. The conservation of stress energy tensor states that 

+ 3-{pi + P,) = 0. (2.2) 
a 

These n+1 equations are sufficient to predict future evolution of the universe, provided 
that the initial value of a and pi are specified. The remaining non-trivial component 
of the Einstein equation gives the dynamical equation 

-2 ^ = 87rG^^P,, 2.3 

but this follows from the previous n + 1 equations. Therefore, it suffices to solve the 
Friedmann equation (12. ip coupled with the conservation equation (12. 2p . However, it 
is also consistent to solve the dynamical equation (12.31) coupled with the conservation 
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equation (\2.2\\ . provided that the Friedmann equation (]2.ip is imposed at an initial 
time. In other words, the Friedmann equation can be considered as an first integral 
of the dynamical equation with a special choice of an integration constant. 

Now, let us suppose that there is a theory without Hamiltonian constraint. Let us, 
however, suppose that in the FRW spacetime, we still have the conservation equation 
(12. 2p and the dynamical equation (12.31) . This is perfectly fine as we have n + 1 
independent differential equations for n + 1 variables, a{t) and pi{t) {i = 1, ■ ■ ■ ,n). 
Actually, we obtain 



as a first integral of the dynamical equation, where C is an integration constant, and 
this is almost the same as the Friedmann equation (12.11) . The only difference is the 
term C/a^. What is interesting is that this is exactly of the form of dark matter. 
In general relativity, dark matter (oc a~^) is included as one of p^'s and, thus, must 
be derived from an action principle since pi is a component of a stress-energy tensor. 
In general relativity this is the origin of the mystery: we need to explain what dark 
matter is made of by specifying its action. On the other hand, in this hypothetical 
theory without Hamiltonian constraint, the term proportional to emerges as an 
integration constant and, thus, we do not need an action for it. 

Intriguingly enough, as we shall briefly explain in the next paragraph, in Hofava- 
Lifshitz gravity there is no Hamiltonian constraint as a local equation at each spatial 
point. Instead, the Hamiltonian constraint equation in Hofava-Lifshitz gravity is an 
equation integrated over a whole space. In homogeneous spacetime such as the FRW 
spacetime, the global Hamiltonian constraint is as good as local one since all spatial 
points are equivalent. However, in inhomogeneous spacetimes there can be drastic 
differences. If the whole universe is much larger than the present Hubble volume 
then it is possible that the universe far beyond the present Hubble horizon is different 
from our patch of the universe inside the horizon. In this case, the global Hamiltonian 
constraint does not restrict the universe inside the horizon. Even if we approximate 
our patch of the universe inside the present horizon by the FRW spacetime, the whole 
universe can include inhomogeneities of super-horizon scales and, thus, the global 
Hamiltonian constraint does not restrict the FRW spacetime which just approximates 
the behavior inside the horizon. Therefore, as in the hypothetical theory considered 
in the previous paragraph, the absence of local Hamiltonian constraint in Hofava- 
Lifshitz gravity results in an extra term oc in the "Friedmann equation" or, to be 
precise, the first integral of the dynamical equation. As before, this term can mimic 
dark matter but we do not need an action for it. 




(2.4) 
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Absence of local Hamiltonian constraint in Hof ava-Lifshitz gravity originates from 
the projectability of the lapse function. The basic quantities in Hofava-Lifshitz gravity 
are the 3-dimensional spatial metric Qij, the shift vector A^* and the lapse function 
A^. In terms of these quantities the 4-dimensional spacetime metric is written in the 
ADM form: 

ds"^ = -N^df + gij{dx' + N'dt){dx^ + N^dt). (2.5) 

The former two, Qij and A^\ can depend on both spatial coordinates x^ and the 
time variable t. On the other hand, the projectability condition states that the lapse 
function A^ should depend only on t and be independent of spatial coordinates. The 
projectability of the lapse function stems from the fundamental symmetry of the 
theory, i.e. invariance under the foliation-preserving diffeomorphism: 

x'^x\x\t), t^i{t), (2.6) 

and therefore must be respected. Essentially, the lapse function A^ represents a gauge 
degree of freedom associated with the space-independent time reparametrization. 
Thus, it is very natural to restrict A^ to be independent of spatial coordinates. This 
is the projectabihty condition. Of course, a space-independent A^ cannot be trans- 
formed to a space-dependent function by foliation-preserving diffeomorphism. This 
point was already made clear by Horava pLj 0- 



3 IR limit of Hofava-Lifshitz gravity 



In the IR limit the action of Hofava-Lifshitz gravity is reduced to 

Ihl = — ^ / dtd^x^N Ik^jK'^ -XR^ + R- 2A] , (3.1) 

where Kij = {gij — DiNj — DjN.j) / {2N) is the extrinsic curvature of the constant time 
hypersurface, K = K], and D is the 3-dimensional covariant derivative compatible 
with Qij. This looks identical to the Einstein- Hilbert action in the ADM form if and 
only if A = L Hence, hereafter, we assume that the renormalization group (RG) flow 
brings A to 1 in the IR or that A stays at 1 from higher energy scales all the way down 
to the IR under the RG flow. The RG flow of Hofava-Lifshitz gravity has not been 
investigated in details and, thus, must be addressed in the future. In this paper, we 
simply assume that A = 1 is an IR flxed point of the RG flow. 



^In the last paragraph of subsection 2.1 of [T], it says that, except for the case with extra symmetry 
such as the Weyl symmetry (A = 1/3), fluctuations of the lapse function must be space- independent. 
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Even with A = 1, however, there is an important difference between the IR hmit 
of Hofava-Lifshitz gravity and general relativity. In Hofava-Lifshitz gravity the pro- 
jectability condition requires that the lapse function should depend only on t. 
Because of this restriction, the Hamiltonian constraint, i.e. the equation derived 
from functional derivative of the total action with respect to the lapse function, is 
not a local equation but an equation integrated over a constant time hypersurface: 

I d'x^iG^il + Kg^} - 87TGNT,,)nV = 0. (3.2) 

Here, gj^J is the 4-dimensional metric shown in fl2.5l) . G^fj is the corresponding 4- 
dimensional Einstein tensor, T^^, is the stress energy tensor, and is the unit normal 
to the constant time hypersurface given by 

n^dx^ = -Ndt, n^^d^ = j^{dt-N'd,). (3.3) 

On the other hand, the momentum constraint and the dynamical equations are local 
equations as in general relativity: 

(gSJ + A^7fJ - 8nGNT,^)n^ = 0, (3.4) 

and 

G\f + Aglf - SttGnT^j = 0. (3.5) 

The Hamiltonian constraint (13. 2p and momentum constraint (13.41) are preserved 
by the dynamical equations ( 13. 5p . Thus, it suffices to solve the dynamical equations, 
provided that the initial condition satisfies the constraint equations. Note that the 
global Hamiltonian constraint (13. 2p is less restrictive than its local version, and allows 
a richer set of solutions than in general relativity. 



4 Dark matter as "integration constant" 



Let us define deviation from general relativity by 

T^u' = ^ (G^it^ + ^9j!J) - T,., (4.1) 

or 

+ Ag^^J = SrcG^ {t,, + Tf/) . (4.2) 

This looks like Einstein equation with the dark sector T^^. In the IR, not only 
the gravitational sector but also the (real) matter sector should respect the full 4- 
dimensional diffeomorphism invariance. Thus, the conservation of energy momentum 
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tensor for real matter W^T^y = should hold in the IR, where V is the 4-dimensional 
covariant derivative compatible with g'^^^. The Bianchi identity then implies the con- 
servation V^Tl^J" = of the dark sector. If the 4-dimensional diffeomorphism invari- 
ance is slightly broken in the (real) matter sector then V^T^^ = — V^T^j^ 7^ 0. Note 
that invariance of the (real) matter action under 3-dimensional spatial diffeomorphism 
implies that V^T^j = and thus V^T^y oc riy. 

The field equations in the IR limit of Hofava-Lifshitz gravity is now written in 
terms of T^^. The Hamiltonian constraint (I3.2p is 

j d'x^Tl^J^n^n' = 0. (4.3) 

The momentum constraint 03.41) and dynamical equations 03.51) are 

T^J^n'^ = 0, (4.4) 

and 

= 0. (4.5) 

As a general solution to the momentum constraint and dynamical equations, we 
obtain 

Tj"/ = P^'^^n., (4.6) 

where p^^ is a scalar function of spacetime coordinates {t,x^). This is equivalent to 
the stress energy tensor of a pressureless dust with energy density p^^ and the unit 
tangent n'^ to its flow. Note that is tangent to a congruence of geodesies: 

n'^V^n, = n^V.n^ = ^d{n>'n^) = 0. (4.7) 

Here, for the first equality, we have used the expression 03. 3p and the fact that the 
lapse function N depends only on t. Finally, the Hamiltonian constraint is 



J d^x^p"^ = 0. (4.^ 



This states that the total energy of the dust-like component in the dark sector should 
vanish. Of course p^^ can be positive everywhere in our patch of the universe. As 
already explained in Sec. [2], this applies even when we approximate our universe by 
a FRW spacetime since overall homogeneity inside the current Hubble volume does 
not exclude super-horizon inhomogeneities. 

As stated at the end of the previous section, the dynamical equations preserve the 
constraint equations. Thus, it suffices to solve the "modified Einstein equation" 

+ Agl!^ = 8nGN {t,. + p^'^n.n,) , (4.9) 
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coupled with field equations of real matter fields, provided that the initial condition 
of the dark sector satisfies the global Hamiltonian constraint (14. 8p . Note that p^^ 
does not have to vanish everywhere. It can be positive somewhere in the universe and 
negative elsewhere, as far as it sums up to zero. For example, p^^ can be positive 
everywhere in our patch of the universe inside the present Hubble horizon. Note 
also that the additional term p^^n^ni, is just an "integration constant" and does not 
represent a real dust. In other words, this additional term acts as cold dark matter 
but does not require an action. 

The Bianchi identity applied to the "modified Einstein equation" (14.91) leads to 

(9^p^^ + TTp^^) = -V^T^„ (4.10) 

where d± = n'^d^ and V is the 4-dimensional covariant derivative compatible with 
g'^^^. This is consistent with the fact that the 3-dimensional spatial diffeomorphism 
invariance of the (real) matter action implies V^T^j^ oc Uy. Therefore, we obtain 

5±p^^ + iTp^^ = n'^V^T^y. (4.11) 

The right hand side acts as a source term for the "dark matter" energy density 
p^^ , and is non- vanishing if the (real) matter sector breaks a part of 4-dimensional 
diffeomorphism invariance. 

If the gravity sector also breaks the 4-dimensional diffeomorphism, i.e. if A slightly 
deviates from 1 or/and higher spatial curvature terms become non-negligible, then 
the left hand side of the "modified Einstein equation" (14. 9 p gets corrections as 

G^^l + Agf^li' + 0(A — 1) + (higher curvature corrections) = SvrGiv {t^i, + p^^n^j^n^^ . 

(4.12) 

The new terms in the left hand side result in extra contributions to the right hand 
side of fl4.10p but they are again proportional to riy due to 3-dimensional spatial 
diffeomorphism invariance. Therefore, in general we have 

d^p"^ + Kp^^ = n^V^T^y + 0(A - 1) + (higher curvature corrections). (4.13) 

The right hand side is non- vanishing and the "dark matter" is inevitably "generated" 
if either the (real) matter sector or the gravity sector breaks a part of 4-dimensional 
diffeomorphism invariance. Thus, turning off the "dark matter" completely or, equiv- 
alently, imposing a local Hamiltonian constraint is not a consistent truncation in 
general. 

In the early universe, the r.h.s. of (I4.13P is non- vanishing. Quantum fluctuations of 
scalar graviton, tensor graviton and (real) matter fields act as source of "dark matter" . 
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Once a cosmological model is specified, it is possible to predict the typical amplitude 
of p^^ generated in the early universe. (Such quantum fiuctuations include modes 
with various wavelengths, including those much longer than the size of the present 
Hubble scale.) For this reason, the initial condition of "dark matter" is not arbitrary 
and this scenario has predictability. 

On the other hand, at late time the r.h.s. of fl4.13p should vanish if matter sector 
recovers 4-dimensional diffeomorphism in the IR and if A = 1 is a stable IR fixed 
point of the RG fiow. In this case, (14.131) is reduced to the usual conservation law for 
"dark matter" : 

d±p"^ + Kp^^ = 0. (4.14) 

The fiow of "dark matter" is tangent to the vector n'^ defined in (13.30 and thus 
is orthogonal to the constant time hypersurface. When a cusp is about to form, 
the spatial curvature of the constant time hypersurface increases. The system enters 
the non-relativistic regime and higher spatial curvature terms become important. 
Among them, terms with z = 3 generate the strongest restoring force. Also, A 
may deviate from 1 by RG fiow. As in some early universe models [H El E], we 
expect that the would-be cusp can bounce at short distance scales. Note that this is 
not because of deviation from geodesies but because of repulsive gravity. In usual 
situation, congruence of geodesies would form caustics and thus cusps because gravity 
is attractive. On the other hand, for the fiow of "dark matter" proposed in the present 
paper, higher curvature terms become important near the would-be cusps and provide 
negative effective energy and repulsive gravity. That is the reason why we expect 
bounce. Note that the bounce is provided by nonlinear terms and thus nonlinear 
analysis is needed. It is also important to include backreactions of the higher spatial 
curvature terms to the geometry since, as mentioned in the previous paragraph, the 
bounce is not due to deviation from geodesies but due to repulsive gravity at short 
distances. Without taking into account nonlinear terms and backreactions to the 
geometry, we would never be able to describe "dark matter" properly. 

The initial value formulation of Hof ava-Lifshitz gravity consists of dynamical equa- 
tion, global Hamiltonian constraint, local momentum constraint and gauge conditions. 
In this language, the "dark matter" emerges only after solving the system of equa- 
tions when we try to interpret a solution. Therefore, as far as scalar graviton, tensor 
graviton and matter fields are properly included as dynamical variables, we do not 
have to consider "dark matter" as an independent dynamical field. This structure 

^ In the case of ghost condensate [15], the derivative of the scalar field responsible for the con- 
densate deviates from geodesies because of higher derivative terms [50]. On the other hand, in 
Hof ava-Lifshitz gravity the vector n'' always satisfies the geodesic equation ()4.7|) . 
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should persist even in quantum level. For this reason, we do not have to promote 
the "dark matter" from non-dynamical integration "constant" to a dynamical field. 
Nonetheless, when we interprete a solution after solving the dynamics of the sys- 
tem either classically or quantum mechanically, we see that difference from general 
relativity amounts to "dark matter" . 

5 Summary and discussion 

In the non-relativistic, power-counting renormalizable theory of gravitation recently 
proposed by Hofava, the so called projectability condition must be respected as it 
stems from the fundamental symmetry of the theory, i.e. the foliation-preserving dif- 
feomorphism. Essentially, the lapse function represents a gauge degree of freedom 
associated with space-independent time reparametrization. Thus, it is very natural to 
restrict to be independent of spatial coordinates. This is the projectability condi- 
tion. Of course, a space-independent cannot be transformed to a space-dependent 
function by foliation-preserving diffeomorphism. The projectability condition then 
implies that the Hamiltonian constraint is not a local equation satisfied at each spa- 
tial point but an equation integrated over a whole space. This point was already 
made clear by Hofava in [1] (See footnote H] of the present paper.). 

Abandoning the projectability condition and imposing a local version of the Hamil- 
tonian constraint would result in phenomenological obstacles [35] and theoretical in- 
consistencies [36j- Note that a strong self-coupling of the scalar graviton reported 
in [35] is not a problem if there is a phenomenon analogous to Vainshtein effect [H] 
since, unlike massive gravity [45J Horava-Lifshitz gravity is supposed to be UV com- 
plete. Other problems reported in [2Sl ES] disappear if the projectability condition 
is respected and if only the global Hamiltonian constraint is imposed. The Poisson 
brackets of constraints form a closed structure since there is only one Hamiltonian 
constraint and it commutes with itself [2] . The divergent coupling of the scalar gravi- 
ton to matter source does not exist in the absence of local Hamiltonian constraint 0. 
In conclusion, both theoretical consistencies and phenomenological viability require 
that the Hamiltonian constraint is not a local equation but an equation integrated 
over a whole space. 

^In the case of massive gravity, Vainshtein effect [44] removes the vDVZ discontinuity [46l but 
simultaneously makes quantum corrections uncontrollable |48j . For this reason, the massive gravity 
theory without a UV completion docs not have predictability even at macroscopic scales. 

^If the local Hamiltonian constraint were not used, then eq. (68) of [35] would not show a divergent 
coupling. This can be seen by moving the first term (written in term of a) in the r.h.s. to the l.h.s. 
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The global Hamiltonian constraint is less restrictive than its local version, and 
allows a richer set of solutions than in general relativity. We have shown that a 
component which behaves like pressureless dust emerges as an "integration constant" 
of dynamical equations and momentum constraint equations. Consequently, classical 
solutions to the infrared limit of Horava-Lifshitz gravity can mimic general relativity 
plus cold dark matter. The "dark matter" satisfies the (non-)conservation equation 
f l4.13p . whose source term is turned on if a part of 4-dimensional diffeomorphism 
invariance is broken by either the (real) matter sector or the gravity sector. In the 
IR, the source term vanishes and the standard conservation equation fl4.14p holds. 
Also, the "modified Einstein equation" (14.121) leads to the Poisson equation (in a 
gauge with N = N{t)). Of course the "dark matter" can cluster. 

Note added 

After the present paper had appeared, the emergence of the "dark matter" was con- 
firmed in several papers [511 [521 [53] . For example, eq. (22) of [53] is exactly equivalent 
to Einstein gravity plus cold dark matter. Obviously, it is not appropriate to inter- 
prete it as a scalar-tensor theory. (This is obvious if we notice that we usually do not 
try to interprete other cold dark matter models as scalar-tensor theories.) Nonethe- 
less, ref. [53] did. Also, the analysis of ref. [53] is not capable of seeing discontinuity 
in the limit A — > 1 since A is set to 1 from the beginning. 

Ref. ^52] has three statements about Hofava-Lifshitz gravity with the projectability 
condition: (i) "dark matter" forms caustics; (ii) "dark matter" is described by ghost 
condensate [IH]; (iii) the scalar sector gets strongly coupled at the scale A ~ p]/^ 
even with A 7^ 1. Actually, these three comments are not correct for the following 
reasons, (i) They did not take into account repulsive gravity due to nonlinear higher 
curvature terms (see the second-to-the-last paragraph in Sec. H] of the present paper 
and also ref. [5^). (ii) Ghost condensate and Hofava-Lifshitz gravity have different 
symmetries. In particular, unlike ghost condensate, the invariance under (space- 
independent) time reparametrization forbids /iqq and thus tt^. (iii) The strong coupling 
away from A = 1 found in [52] indicates breakdown of their description, i.e. the way 
dynamical degrees of freedom are identified, but does not imply inconsistency of the 
underlining UV theory. Note that the scalar graviton does not get strongly coupled 
away from A = 1 as is clear from e.g. eq. (54) of ref. [T]. See also the last paragraph 
in Sec. IHof the present paper. When we take the limit A — 1, we have to take into 
account nonlinear interactions carefully to see if there is an analogue of the Vainshtein 
effect t44j. 
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More recently, ref. ^55] rediscovered ghost instability in the regime 1/3 < A < 1. 
This can be seen already in eq. (54) of [T] and just implies that we need to consider 
the region A > 1. 
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